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Abstract. Wc present a residual-based a posteriori error estimate for the 
Electric Field Integral Equation (EFIE) on a bounded polyhedron Q with 
boundary F. The EFIE is a variational equation formulated in H^^^ (F). We 
express the estimate in terms of L^-computable quantities and derive global 
lower and upper bounds (up to oscillation terms). 

1. Introduction 

The Electric Field Integral Equation (EFIE) describes the scattering of electro- 
magnetic waves on a perfectly conducting obstacle fl with surface T, in our case a 
polyhedron. Assuming a time-harmonic dependence, the Stratton-Chu representa- 
tion formula expresses the electric field E in terms of a surface potential as 

E{x) = E'"%x) + ^ (^G,(x,y)u{y) + ^gradr^,G,(a^, y)divrw(y)^ da{y), 

where k denotes the wave-number and is the given incident wave that is 

scattered on T. Invoking the boundary condition that the tangential component 

of the total electric field E vanishes on the surface T, as corresponds to 11 being 

perfectly conducting, the EFIE consists of seeking the surface current u € iJ^iv (^) 

such that for all x E T 

(Gk{x,y)u{y) + pgradr,^Gfe(a;, y)divrM(y)^ da{y) = -7||(£;"'=)(a;), 

where H^^^ (F) is the space of traces of ff(curl , ft) functions that are rotated by 
a right angle on the surface and 7 y denotes the tangential trace onto T. 

Computing approximations of the EFIE by means of the Boundary Element 
Method (BEM) , namely using a Galcrkin approach based on the variational formu- 
lation of the EFIE, is expensive due to the dense matrix structure of the ensuing 
linear system. Although fast techniques such as the Fast Multipole Method exist, 
c.f. [15] as an example of a first work in this field, it is still crucial to locate the 

degrees of freedom efficiently, namely in regions of low regularity of the solution u. 

_ 1 

Since u S ^^^jiv i^)' ^ exhibits in general rather low regularity and, as a conse- 
quence, a priori estimates show extremely low convergence rates for quasi-uniform 
mesh refinements; see [10, 16]. In contrast, adaptive refinement techniques, based 

R. H. NocHETTO, Department op Mathematics, University op Maryland, College Park 
MD 20742 (rhnSmath. umd.edu) 

B. Stamm, Department op Mathematics, University op Calipornia, Berkeley and 
Mathematics Department, Lawrence Berkeley National Laboratory, Berkeley 94720 
(stammSmath . berkeley . edu) 



2 A POSTERIORI ERROR ESTIMATES FOR EFIE ON POLYHEDRA 

on a posteriori error estimates, exploit much weaker regularity of it in a nonlinear 
Sobolev scale and allow for optimal error decay in terms of degrees of freedom in 
situations where quasi-uniform meshes are suboptimal. The design and analysis of 
a posteriori error estimators is, however, problem dependent; we refer to [9, 17] for 
an account of the theory of adaptive finite element methods in the energy norm for 
linear second order elliptic partial differential equations in polyhedra. 

The a posteriori error analysis and corresponding theory of adaptive mesh re- 
finements for BEM is much less developed and an overview of different approaches 
for the former is given in [7, 13, 18]. It seems that this is the first contribution with 
specific application to electromagnetic scattering problems on polyhedra. 

For integral equations, additional difhcniltics arise since the residual typically 
lies in a Sobolev space with fractional index that is possibly also negative, as in 
the present case. Since such norms are not computable in practice, this imposes 
additional challenges to the residual based approach of a posteriori error estimates. 

In this paper we develop nevertheless a residual based a posteriori error estima- 
tor for the EFIE on polyhedra, and prove upper and lower global bounds. Residual 
based estimators are especially attractive due to their simplicity of derivation and 
computation, but they involve interpolation constants which can at best be esti- 
mated. Alternative estimators have been proposed, mostly for elliptic problems 
defined in il, at the expense of their simplicity; we believe that our approach can 
be extended to those estimators as well. We derive computable L^-integrable quan- 

titles to estimate the error of the BEM measured in the H^^^ (F) norm, which is the 
natural norm for EFIE. We therefore avoid evaluating fractional Sobolev norms. 

For proving well-posedncss of the exact solution and developing a priori error 
estimates it is important to decompose both the exact solution and test function 
using a Helmholtz decomposition [6, 16]. In contrast, to derive a posteriori er- 
ror estimates, it is crucial to decompose the test function according to a regular 
decomposition which extends the Helmholtz decomposition; see [8] for i?(div;f2). 

This paper is organized as follows. In Section 2 we recall the necessary functional 
analysis in order to derive a posteriori error estimates for the EFIE [2, 3, 4, 5]. 
We also present and study a Clement type interpolation operator for the Raviart- 
Thomas space, based on ideas from [1]. We discuss the EFIE integral equation in 
Section 3, and derive global upper and lower a posteriori error estimates in Section 
4. Section 5 is finally left for conclusions. 

2. Functional Spaces and Differential Operators 

The functional analysis framework developed in [2, 3] will be used in this work. 
In this section we give a short introduction to the functional spaces and differential 
operators used in the following sections. However, for a detailed and thorough 
overview we refer to [2, 3, 5, 6]. References [5, 6] deal with non-smooth Lipschitz 
surfaces, thus the theory is also valid for polyhedra, and covers therefore a more 
general framework. However, we restrict our theory to polyhedral surfaces. 

2.1. Spaces, norms and trace operators. Let be a bounded polyhedron in 
M^, and denote its boundary by F and its different faces by Fj, j = 1, . . . , Np. The 

exterior part $7+ is defined by r2+ = M'^\0. Let n(x), x E T, denote the outer 
unit normal to the surface F, which is piecewise constant on F. We also indicate 
by Bij = dVi n dFj the edges of F and by the unit vectors parallel to e jj , with 
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its orientation fixed but arbitrary. If = nlr^, we furtlier define 

Ti — X Tti^ Tj — Tij X Tlj 

to be unit vectors lying on the supporting planes of Fj and Fj ; see Figure 1 for an 
illustration. 




Figure 1 . Local coordinate systems around an edge Cij = dTi n dVj . 

On F, we define the space of square intcgrablc tangential fields 

Ll{T) = {v&[L^{V)f\v-n = a.e.}. 

Moreover, we let H\V) and H^{r) = [H^{T)Y, with s G [-1, 1], denote the stan- 
dard Sobolev spaces of complex-valued scalar and vector-valued functions on F and 
denote their norms by || • ||H»(r) and || • ||jfs(r), respectively; for negative Sobolev 
indices the norms are defined by duality. Furthermore, for s G (0, 1), we denote by 

7 : H'+'^ (O) -^H\V), 7 : [H'+ s H' (F) 

the standard continuous trace operators, and by and Ity their continuous right 
inverses. 

For complex- valued vector functions we introduce the facewise 5 -broken space 

m{T) = {v&Li{T)\v\r,&H"^Ti),l<i<Np}. 
with corresponding norm 

|2 i|„.||2 



Ell- 

Moreover, we define the spaces 

- f - 1/2 1 

^\\ (r) = N e -H"- (F) v\ri ■ Tij = v\rj ■ nj, for every edge Cij V , 



^ r ^ 1/2 1 

^K^) = ^ HI (F) v\Yi ■ n = v\rj ■ Tj, for every edge \ , 

where the relation 'I' is understood in the sense that 

(2) v,^l\, ^ j j ^^^M^^d.(.)d.(y)<oo. 

JVi 11*^ y\\ 
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We further define 

■■= ||a._^||3 da(y), 

ir. ip. ^ ||a;-y||3 da(a;)da(y), 

for each edge e^j of the polyhedron and denote by Ij the set of indices i such that 

Tj and have a common edge eij. 

Proposition 2.1 ([2, Prop. 2.6]). The spaces H^^ (F) and H]_{r) are Hilhert spaces 
when endowed with the norms 

Nf 

(3) \\vf r ■.= \\vf r +EE-^'lW' 

In other words, v G if (F) . (F) satisfies v G if 2(rj) on the faces F^ of F, 
and the paraUel resp. orthogonal component of the function v to edges Cj^ of F are 
"if 2 -continuous" in the sense of (2); c.f. [2] for further details. 

We denote by H~^r), H~^r) the dual spaces of -ff|^(r), ifl(r) with pivot 
space (F). The corresponding duality pairing is denoted by (•, •)||^r resp. (•, •)_l,f. 

The norms |j • || _i and || ■ || _i are defined by duality. 
^ (p) ffj_ ^ (p) 

For complex- valued functions v G [C°°{Q)]^ the tangential traces are defined by 

(5) 7||(t;) := n X X n)|j,, 'yj_{v) := (v x n)\^. 

We point out that 7 y (v) = v — {v ■ n)n gives the component of v tangential to F, 
whereas 'y±(v) provides a tangent vector field perpendicular to 7||(t)). Since F is 
a polyhedron, for any edge e,j of F the components of 711(1?) and 7_l(w) tangential 
and normal to e^j are continuous, namely, 

(6) 7||('w)|Pi • Tij - 7 II W If,- • nj = 0, 'y±{v)\r, ■ n - 7-LWIr^ • tj = 0. 

This means that both operators 7 || and 'yj_ can be viewed as face-by-face projec- 
tions; see [2, p. 36]. Combining this observation with definitions (1), we realize that 

i?| (F) and H^{T) are the trace spaces of 7||, 7_l acting on H^{Q). This is stated 
in the following Proposition. 

Proposition 2.2 ([2, Prop. 2.7]). The trace operators 

7 II : H\^) ^ ifjf (F), 7^ : H\n) ^ ff|(F) 

are linear, surjective and continuous operators. In addition, there exists continuous 
right inverse maps : (F) i?i(0) and R± : Hl{T) ifi(Q). 

We can now establish a critical result for the upcoming analysis. Note that 
1 1 
iy^(r) consists of tangential vector fields whereas H^(r) does not. 
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Lemma 2.3. There exists a continuous map t± : H^iT^ — >• H]_(T) with right 
inverse tj^ : i?|(r) ^ Hi{r). 

Proof. We define t±:Hi{r)^ Hl{r) and t^^ : Hl{T) Hi{T) by 

t±{w)=^±{iiy{w)), VweJf5(r), 

where 7 and are the trace and its right inverse, whereas 7_l and -R_l are the op- 
erators of Proposition 2.2. The continuity of these operators implies the continuity 
of t± and tj^- 

To prove that tj^ is the right inverse of t±, we observe that 

t±{tl\v)) = 7±(i^(7(i2±W))), Vt; e HliT), 

and that 7j^ projects face by face on F [2, page 36]. li w = tj^{v) E Hi{T) and 
g = RyW e H^{Q), then w = 'y{R±{v)) and t±{w) = 7_L(g). Since 'y±{g)\ri = 
7(fif)|r ^ for sach face Fj of F, we obtain on Fj 

t±{tj^{v)) = t±{w) = 7_L(t/) = 7(9) X n 

= ^{R^w) X n = w X n = 'y{R±{v)) x n = 'y±{R±{v)) = v. 

Thus, tj^ is indeed the right inverse of tj,. □ 
2.2. Tangential differential operators. We set -^'(F) := 7(iJ^(f2)), and define 

3 1 si 

the tangential operators gradp : if 2(F) H^^{r) and curlr : -ff2(F) i?|(F) 

by 

(7) gradp(^ := 7 y (grad 0) curlr<?f> := 7_L(grad 0) e //"^(fi), 

where grad denotes the standard gradient in M^. According to definitions (5), 
gradp(/) is the orthogonal projection of gradcf) on each face of F, whereas curlr^ 
is obtained from the former by a tt/2 rotation. It can be shown that the maps 
gradp and curlp are linear and continuous. 

The adjoint operators divp : H~^T) H-i{T) and curlp : H~Ht) 

H~ 5 (F) can be defined as follows 

(divrv,(^)|_r = -(u, gradp?!') ||,r, 
(curlptu,(;i)| P = {w, curlr (t))±,r, 

for all (j) e Hi{r), v G -H'P(F) and w e H~Ht). 

In view of these definitions we now introduce the spaces 



-^curl 



J(F) := G ^(F) I divrv G ^^-^(r)} , 
i{T) := {v G H--^r) I curlpt; G ^^-'(F)} , 
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which wiU play a crucial role for the upcoming analysis and are endowed with the 
graph norms 

\\vf , := \\vf 1 +||divrt;||' i, . 
"ff;j(r) ^F^^(^) "ff-^r)' 

\\v\\^ 1 := 11-^11' 1 +||curlr'y|P i . 

Let the natural solution space of Maxwell's equations be denoted by 
i?(curl ,n) = {v e L^{n) I curl v e L^ifl)} . 

Theorem 2.4 ([3, Theorem 4.6]). The mappings 7|| and 7_l admit linear and 
continuous extensions 

7 II : if (curl , Q) ^ Jf;i(r), 7± : (curl , Q) ^ H^^ (F). 

Moreover, the following integration by parts formula holds true: 

J ^curl t; • M — curl M • t;^ dn = (7_lW, 7||t;)||^r, ^ u G H {curl ,fl),v G H^{fl). 

_ 1 _i 

Furthermore, a duality pairingj_(-, •)|| between H^^^ (F) and H^^^i{r) can be es- 
tablished by using an orthogonal decomposition of those spaces so that the following 
integration by parts formula still holds 

J ^curl u • u — curl w • dfi =_l(7_lu, 7111;) II , y u,v G H{curl 

For more details, we refer to [3]. 

The differential operators gradp and curlp can be further extended as follows. 

Proposition 2.5 ([3, p. 39]). The tangential gradient and curl operators introduced 
in (7) can be extended to linear and continuous operators defined on H2{r) 

gradr : H^T) ^ H~^T), curlr : h'^T) ^ H~^T). 
Moreover their formal i|(F)-adjoints 

divr : (F) ^ H' 3 (F) and curlr : (F) ^ H' s (F) 
can be defined by 

(divrv,(/>)i r = -{v,grad^(l))±,r, 
(curlrii;,(/))i_r = {w, curlr (t)}\\^r, 

for all (pG Hi{T),v G HI (F) and w G h'^ (F). 

2.3. Potentials. Let Gk denote the fundamental solution of the Helmholtz opera- 
tor A + A;^, which is given by 

,^ exp(^fc|.-^|) _ 

The scalar and vector single layer potential are then defined respectively by 
(F) ^ HU^% ^^tivKx) := j^Gk{x,y)v{y)da{y). 
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These potentials are known to be continuous [6, Theorem 3.8]. Finally the scalar 
and vector single layer boundary operators are defined by 

Vk:H-^T)^Hi{r), 14:= 70*^, 

The simultaneous continuity of the trace operators 7, 7 y and the single layer po- 
tentials yield then the continuity of the single layer boundary operators, namely, 

(9) \\Vkv\\ i,,^\\v\\ 1,,, \\Akv\\ 1 ^ llull 1 , 
II Ii/f2(r)-" "H-2(r)' " " "ff2(p) - II "ff-2(r)' 

for all V G H^2{T) and v G H^^ ^(r). In particular, if restricted to L'^{T) C 
H-i{r), then the range of Vk lies in H^F) [6, Theorem 3.8], i.e., 

(10) ImiVkiLHT)) C HHT). 

Likewise, for the vector case the corresponding result reads [4, Prop. 2] 

(11) im(Afe(i2(r)) c H\r). 

2.4. Interpolation of weighted spaces. In the following section we will be con- 
fronted with interpolation of weighted i^-spaces. We thus recall in this section 
some basic results taken from Tartar's book [20], which are valid without regularity 
on the weights. 

Let T be a family of shape-regular triangulations decomposing F into flat trian- 
gles such that the surface covered by the triangles coincides with F. Denote the set 
of edges of the mesh by S-r- For a fixed triangulation let hr denote the diameter of 
any element T gT and let h be the piecewise constant function such that h\T = fiT- 

Lemma 2.6 ([20, Lemma 22.3, p.llO]). If A is linear from Eq+Ei into Fq+Fi and 

maps Eq into Fq with \\Ax\\Fa < Mo||a;||Ej, for all x Eq, and maps Ei into Fi with 
\\-^3:\\fi < -Milja^llBi for all x € Ei, then A is linear continuous from {Eq,Ei)0^p 
into (Fq, Fi)g^p for all 9,p, and for < 6 <1 one has 

||Aa||(Fo,Fi).„ < ^^o"'<l|a||(£o,Bi)..p for all a G {Eo,Ei)e,p. 

The space {Eo,Ei)g^p denotes the interpolation space between Eq and E\ based on 
the LP -inner product. 

Lemma 2.7 ([20, Lemma 23.1, p. 115]). For a measurable positive function w on 
T, let 

E{w) = ^ J \u{x)\'^w{x) dx < 00^ with \\u\\w = (^J \u{x)\^w{x) dx^ . 

If wq,wi are two such functions, then for < 9 < 1 one has 



{E{wo), E{wi))g^2 — E(wg) with equivalent norms, where wg = Wq Wi- 

Corollary 2.8. Let < s < 1 be arbitrary. Let A be a linear continuous map from 
L2(F) into L2(r) and from H^{T) into L'^{T) with 

||A^;||i2(r) < Mo||v||i,2(r) for all v e L^(F), 
||h~^^i)||i,2(r) < Mi||ti||Hi(r) for all v e -H'^(F). 
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Then A is a linear map from H%r) = {H^{T), L'^(r))s.2 into L'^{r) with 

||h-M«|U2(r) < M^-'MIMhs^d for all v € ff^(r). 

Proof. Combine Lemma 2.6 with Lemma 2.7. □ 

2.5. Discrete spaces and interpolation. Let RTq{T) denote the /oca^ Raviart- 
Thomas space of complex- valued functions on T e T defined by (cf. [14, 19]) 

RTq{T) := {v{x) = a + /3x I a e C^ ^ e C} . 

The global Raviart-Thomas space is defined by 

RTo ■■= {v e Hl^{T) I v\T e -Rro(T) VT e r} , 

where ilf^;^(r) is defined in a standard manner 

Hl^{T) := [v e Ll{T) \ divrt; e ^^(r)} . 

Further denote by Vr the space of scalar complex-valued continuous functions that 
are piecewise linear, namely 

(12) v(r) = {«€ i/i(r)|«|K €Pi(T)}, 

where Pi (T) denotes the space of affine polynomials on T. Let M{T) denote the 
set of all nodes y oiT and {'^v}u^{t) be the family of nodal bases of V(T). 

Definition 2.9. Let the Clement type interpolation operator Iq- : L'^iT) — )• V(T) he 

lrv~ Mv)<fi^ yv€L\r) 

where = supp((^,^) and the degrees of freedom are given by 

Proposition 2.10. If v £ Jf*(r) with < s < 1, then the interpolation operator 
Ij- satisfies the following interpolation properties 

(13) \\h-'{v - Ir^^)|U2(r) ^ Mn^iD for all v e H'{T). 

Proof. This interpolation operator Ij- is also used in [12] and the following result 
is proven 

||h"i(z; - (/i„(i;))||i2(r^) ^ ||gradrw||i,2(r„) 

for V G H^{T) under the assumption of shape-regularity [12, (2.2.29)] (Note that 
in our c;ase the mesh matches the surface and therefore equation (2.2.29) can be 
simplified). Following the arguments of the original paper of Clement [11, Proof of 
Theorem 1], it is now straightforward to prove that 

\K\v -lrv)\\L^T) IISJ^adrw||i,2(r^), 

where J\f(T) denotes the set of nodes of the element T. Now, summing over all 
elements of the mesh T and using that the number of elements sharing a node is 
bounded, as a consequence of shape regularity of T, we get 

||h"^(t;-Iri))||L2(r) ^ \v\m{r)- 
Furthermore, the operator can also be shown to be i^-stable [12, (2.2.33)]. 
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Figure 2. Illustration of the normals on an element Tg. 

Therefore, the linear continuous operator Aj- = Id — 17- : i^(r) 1-^ -^^(r) satisfies 
||h"^('!;-Irt')||i,2(r) ^ llfllHHr), 

\\v-lTv\\mr) ^ \\v\\mr)- 
The asserted estimate (13) follows from Corollary 2.8. □ 

Besides this for s = 1/2, we will also need a Raviart-Thomas type interpolation 

1 

operator for functions in v G H^{T). Since divrv ^ -^^^(r), the standard degrees 
of freedom are no longer well-defined. Therefore, we will utilize an interpolation 
operator similar to that introduced in [1] for the first type Ncdolcc elements. 

For any edge e e f 7- of the mesh we associate an arbitrary but fixed element Tg 
such that e C dTe- On Tg we denote by TTg the L^(Te)-projection onto constant 
functions. We let {'0e}ee£r ^'^ the standard Raviart-Thomas basis of lowest order, 
sometimes also referred to as the Rao- Wilton- Glisson (RWG) basis in this context 
of electromagnetic scattering, such that 

■j/je • i^e ds = 1 and J xl)^' ■ i^e ds = 

for any e' S £j- such that e ^ e' and where Ve denotes the outer unit normal of Tg 
at the edge e which is coplanar with T^; see Fig 2.5. 

Definition 2.11. Let the Clement type interpolation operator Ij- : Lf{r) RTq 
for the Raviart-Thomas element of lowest order be given by 

eeSr 

where the degrees of freedom are defined by 

ae{v) := j TTeiv) ■ ds. 

Remark 2.12. Note that I7- does not satisfy the usual commutative property 

divr(Ir^^) 7^ Po(divrf ), 

where Pq denotes the element-wise -projection of degree 0. This is important in 
the a priori analysis but not in the upcoming a posteriori error analysis. 

Lemma 2.13. The degrees of freedom of the interpolation operator I7- : L^^T) — > 
RTq are well-defined and I7- satisfies the local L"^ -stability bound 

Utv\\l^{t) < c ||w||i,2(Ar) for all T eT, 
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where At denotes the set of elements that share at least one edge with, T. 

Proof. The argument is similar to [1] . If T e T is an arbitrary but fixed element 
and £{T) denotes the three edges of T, then 

||Ir'w||l-2(T) = II 5Z 0'eV'e||l,2(T) < ll"eV'e||i2(T) = ^ jofej ||l/'e||l,2(T)- 

ei££(T) eeSiT) eeS{T) 

Invoking the Piola transformation, it can be shown that 



HeU^T) < C||l/'e|li,2(f) 



< c 



since the basis functions ^pe on the reference clement T are uniformly bounded 
independent on T. Moreover, applying the Cauchy-Schwarz inequality, we get 



J e 



< c/l|, ||7reV||l,2(e) < chT\\-keV\\L'^(e) 



where TTg denotes the L^(T)-projection onto constant functions on the reference 
element T. Note that TTei) is defined via the affine transformation from Tg (and not 
T) to T. Norm equivalence of polynomials (constant functions in this case), the 
L^-stability of TTe and a scaling argument yield 

(14) |ae| < C/IT II ii-e* II i2(f) < c/lT||t'||i,2(f) < c||^;||i,2(TJ. 

Combining the above estimates implies the asserted stability bound of I7-. □ 

To explore the accuracy of the interpolant I7-, we need the following lemmas. 
Lemma 2.14 (Local approximability). For any c G M.^, we have the error estimate 

\\v - ItvWlht) ^ 11^- - 1J-c\\l^At)i for all T eT, 
where At denotes the set of elements that share at least one edge with T. 

Proof. Let v = jj_c with c e which, in view of (5), is piecewise constant in T. 
According to (6) the normal component of v is continuous across all edges of the 
mesh, including those of the polyhedron F, whence v € RTq. 

We first observe that I-fV = v because tTsvIt^ = v\t^ for any edge e C dT and 



ae{v) = J TTeV • I>'e ds = J V ■ ds. 



Since these three local degrees of freedom on T e T are unisolvcnt and they coincide 
for both lrv\T and v\t, we deduce ItMt = v\t- Consequently 

\\V - ItvWl^T) < \\V - v\\l^T) + Uriv - v)||i2(T). 

By the local L^-stability of Lemma 2.13 we conclude that 

\\v - Itv\\l^{t) ^ \\V - v\\L^{AT)y 

as asserted. □ 

1 1 
Lemma 2.15. If v G H^{T), then there exists w G 2(r) such that t±{w) = v 

and \\w\\ 1 , ^ llvll 1 

Proof. Simply set w = tj^{v), where tj_^ is defined in Lemma 2.3, and use the 
facts that tj^ is the right inverse of t± and the continuity of tj^- 1^ 
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Lemma 2.16. If w e if5(r), then !|t_L(i«)||i,2(AT) ^ \\'w\\l'^{At)- 

Proof. As in the proof of Lemma 2.3, let g G i?^(f2) be the function such that 
g = FLyW and tj_{w) = j^ig). Therefore 

7_l(9)|t = 7(9) It X n, for a.e. x G T, for all T gT, 

whence 

ll*±Wlli2(A^) = J2 \\t^M\\lHT)= \h(9) X n\\l2^T) 

TcAt TcAt 

^ E ll7(9)lli2(T) = ll7(9)lli.(A.) = lkllL(A.) 

TcAt 

because FLy is the right inverse of 7 and thus 7(g) = w. □ 

1 

Proposition 2.17 (Global approximability) . Ifv G H^{T), then the interpolation 

operator I7- satisfies the following global error estimate 

(15) \KHv-'tTv)\\mr)^M 1 forallv&HliT). 

■"_L I 

Proof. Given v G H^{T), there exists w G i?2(r) so that t_L(iu) = v according 

to Lemma 2.15. For each T G T, we define wt = w{x)dx G and vt = 

1 

t_\_{wT) G H^{r). Since the estimate of Lemma 2.14 is local, we have 

||h"(t;-Irt.)||i.(r) = ||h"(t.-Irf)||i.(T) < E W^riv -vt)\\1.^^,) 

for a = —1,0. Using Lemma 2.16 yields 

\\v - VtWl^At) = ~ «'t)||i,2(At) ^ - ^«T||l,2(Ar) 

and stability of the L^-projection together with the definition of Wt implies 

\W -Wt\l^{At) ^ \\w\\lhAt), 
\\h:^\w-WT)\\LHAT) ^ Ikllffi(AT). 

whence 

||u - Irt^||i2(r) < E ll'^lli^(AT) ^ ll"'lli=(r). 

tgt 

||h-i(t;-Irt;)||i2(r) ^ E IMh^a^) ^ IMh^ty 
TeT 

Applying Corollary 2.8 to vector-valued functions, we obtain 

where the last inequality results from Lemma 2.15. This concludes the proof. □ 
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3. Problem Setting 

The variational formulation of the Electric Field Integral Equation (EFIE), also 

_ 1 

called Rumsey principle, consists of seeking u G H^-^ (F) such that 

(16) aiu, v) =i(/, v) II for all v e (F) 
where / € H^^^^(T), the sesquilinear form a(-, •) is given by 

a{u,v) := (\4divrw,divr'u)i^r - k'^{AkU,v)\\^r, 
_ 1 _ 1 

•)|| is the duality pairing between H^^^^{T) and H^^^ (F), (•, •) i^r the duality 

pairing (T) - H'i (T), (•,-)||,r is the duality pairing H^{T)- H~^T), and the 

integral operators Vk-Ak has been defined in §2.3. 

The discrete formulation reads: find U € HTq such that 

(17) a{U, V) =±{f, F>|| for all V e RTq. 

Equation (16) is well-posed under the assumption that the wave number k does 
not correspond to an interior eigenmode of the Maxwell problem on F. As a con- 
sequence, the following continuous inf-sup condition holds (see also [16]): 

(18) _i -< sup „ "^"'""^ forallwei?"|(r). 

^ ' " ' H..^r) ~ 1 \\v\\ 1 divV ; 

Since the boundary element discretization is conforming, i.e. RTq C H^^^ (F), the 

_ 1 

following Galerkin orthogonality holds: if ti e i^^iv (^) solution of (16) and 

U e RTq is the solution of (17), then 

(19) a{u -U,V) = Q for all V e RTq. 

In addition, as a direct consequence of the Cauchy-Schwarz inequality and the 
continuity (9) of the single layer boundary operators, the form a(-, •) is continuous: 

(20) a{v,w) ^ \\v\\ . \\w\\ . for all v,w G H-^{V). 

■"div I -"div (1 ) 

4. A Posteriori Error Analysis 

As is customary in the theory of a posteriori error estimation, one has to assume 
a higher regularity of the right-hand side than it is needed for well-posedness in 
order to derive computable error bounds. Therefore we assume in this section that 

(21) /e/f|f(F)n/rLi(r) 

with (F) given in Proposition 2.1 and if°uri(r) = [v & ^^(F) | curlr?; S i?(F)}. 

We proceed as in Gascon, Nochetto, and Siebert [8] for flat domains. To this 
end, we start with some auxiliary results that will be useful for our analysis later. 

Lemma 4.1 (Regular decomposition [5, Theorem 5.5]). The decomposition 
H-^ (F) = curlr(F^ (r)/C) + ffl(F) 
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is valid and is stable, i.e. for v — + curlrQ; with G H]_ (F) and a £ (r)\C, 

(22) 11*11^^^^^^ + ||a||^.^^^ < ll^ll«-^i(r) ^ ^ ^^''^ 
Lemma 4.2. ForY{T) given by (12) there holds curlr(V(T)) C RTq. 
Proof. By [5, Corollary 5.3] we have 

ker(divr) n Z?(r) = curlr(F^(r)). 

Thus for all a S V(T) C H^(r) wc infer that curlpa € Lf{r) is piecewise constant 
and that divpcurlra = G i^(r). This implies that curlra G iJ°i^(r). □ 

_ i 

4.1. Upper Bound. Let u G H^.^ (F) be the exact solution of (16) and U G -RTq 
be its approximation defined by (17). By the Galerkin orthogonality observe that 

a{u -U,v)= a{u -U,v-V) for all v G iJ^il (r), V G RTq. 

Decompose v as v = ^ + curlra, according to Lemma 4.1, and define 

(5* := * - *r, 6a:=a-ar 

where '^'7- G iJTo and ctT- G V(T^ can be arbitrarily chosen. Thus we can write 
V — V = + curlpfc and 

a{u -U,v-V) =±{f, + curlrfo) || - a{U, (5* + curlpfo) 

=±(/,'5*)|| +(fc'AftC/,'5*)||,r 
^ </ ' 

+_L(/,curlrfo)|| + (A:^AfcC7,curlr(5a)||,r - (VfcdivrC/, divr(5*) i r 
^ . ' V ^ '-^ 

=I2 =13 

_ 1 

for any v G -ff^iv (r)' *T ^ -^^o and aj- G V(T). We proceed in four steps. 

Wc note that / G i.?(r), fc^^feLf G i?|f (L) C Lj{T) and that * G i?l(r) C 

i|(r) due to enhanced regularity of * asserted in Lemma 4.1. Since *r S -R^o C 
(r), we can replace the duality pairing in Ii by an integral and thus write 

(23) Ii= J {f + k^AkU)-6^da. 

Since / G H^^ (F) the duality pairing_L(-, •)|| can be interpreted as 

_L(/,curlr(5a)|| = (/, curlrfo)||,r, 

1 

namely as a duality pairing in H^^ (F). The definition (8) of curlr now yields 

Ii = {f + k'AkU, curlr<5a)||,r = (curlr(/ + k'AkU), 5a) 

Since 5a G Hh{T) and curlr(/ + PAkU) G ^^(r) because of (11) and (21), we 
can also write I2 as an integral 

(24) "^^"Z curlr(/ + A;^AfeC7)&5;dCT. 

[3] For the last term I3, we integrate by parts according to (8), whence 
I3 = -(gradr(yfedivrC/),5*)^,r. 
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Since divr(-RTo) C ^^(r) we infer that gradr(Vfcdivrf7) e L?(r) in light of (10). 
This impUes that I3 is also an integral 

(25) = - y" gradr(Vfcdivri7) • 5* dcr. 

Inserting (23)-(25) back into the sesquilinear form a yields 

(26) a{u-U,v) = j R-d^Aa + j rSada for all u e -Ff^| (F), 



IT Jr 



where R G Lf{T) and r e L (T) are given element-by-element by 

R\j, := f + k'^AkU + gradr(VfcdivrC/) for aU T e T, 
^^^^ r\T ■■= curlr(/ + k^AkU) for all T e T. 



We now choose aj- = l-j-a and ^^7- = where I7- and I7- are the interpolation 
operators of Definitions 2.9 and 2.11. Applying the Cauchy-Schwarz inequality and 
the interpolation estimates (13) and (15) yields 

a{u-U,v) < ||h5i2||i2(r)||h"3(5*||i2(r) + ||li3r||i2(r) ||h-3fo||i,2(r) 
^ lMH|U2(r)||*||^|^^^ + ||h*r|U2(r)||a||^i^^^, 

which together with the stability (22) of the regular decomposition leads to 
a{u-U,v) ^ (||h5l?j|i2(r) + ||h5r||L2(r)) \\v\\ 



Combining this with the inf-sup condition (18) finally implies 

ki-J7 _i -< sup -J—, ^ ^ hsil r + li^r- r- 

" "f/,.2(r) ~ 1 \\v\\ _i ~ 

We summarize this derivation in the following theorem. 

Theorem 4.3 (Upper bound). Let f e {T) nH°^^i{T), u e H~^{T) be the 
exact solution of (16) and U € RTq be its approximation defined by (17). Then, 
there exists a constant Ci > depending on shape regularity of T such that the 
following bound holds 

where the element indicators r]q-{T) are defined as follows in terms of the residuals 
R e if (r) and r e L'^{T) given in (27) 



rfriT) ■.= hT\\R\\l.,r)+hT\\r\\^ 



Remark 4.4 (Trace regularity of an incident plane wave). It the right hand side / 
is the tangential trace of a plane wave Einc then we conclude from the analyticity 
of the plane wave and of all its derivatives that 

/ = 7 II (£;,„,) e ff|f (r), 7 y {d^.Ei^,) e ff|f (r) 

for i = 1, 2, 3. Therefore, / satisfies the stated regularity assumption (21). 
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4.2. Lower Bound. We next show a global lower bounds for the error indicators 
r)^{T). Since R € L^{T) and r € i^(r) we define the local constants 

Rt= R{x) da{x) rr ^ r{x) da{x), for all T e T, 
Jt Jt 

and their global piccewise constant counterparts Ro\t — Rt and tqIt = tt- 

_ \ 

Theorem 4.5 (Global lower bound for the residual). Let u G -ff^iv (^) exact 
solution of (16) and U G RTq be its approximation defined by (17). Then, there 
exists a constant C2 > 0, only depending on shape regularity of T, such that the 
following bound holds 

C^W^^RU^iv) < \\u-U\\ , +\\hHR-Ro)\\LHT). 

■"div ) 

Proof. Let 6t : O ^ K be a bubble function, namely a Lipschitz function so that 

supp 6t C T, / bTdx = \T\p:i / b'^dx, 
Jt Jt 

for a given T G T. Such a function can be given by a polynomial of degree three on 
T consisting of the product of all three barycentric coordinates times a real scaling 
factor. Let ^t = CTbT with ctt € and note that as a direct consequence of the 
first point we have 

(28) / divr*T da = vE't • nr ds = 0. 

Jt JdT 

For the particular choice <tt = hTRT, we see that 



i 



Rt ■ ^Tdx = hT\\RT\\L^(T) 

JT 

and 

||*t||i,2(T) ^ hT\\RT\\L^{T) ^ ||*t||z,2(T)- 

We construct a global function ^ so that its restriction to T coincides with ^t for 

1 

all T € T- We claim that * € H]_ (F) because it is made of piecewise polynomials 
with vanishing normal component on the interelement boundaries of T. In view of 
Lemma 4.1, such a * is an admisible test function in (26) and, together with the 
choices ^7- = and a = aj- = 0, yields 

a(M - C/, *) = / H • * do- = / (il - ilo) • * dcr + / Rq ■ "9 da 



{R-Ro) ■ *d£7+ ||h5i?o||i2( 



(D- 



By the continuity of the sesquilinear form a(-, •), we have 
||h5fio||i2(r) = a{u - J7, *) - J^{R - Jio) • * dx 



(29) 



11*11 -i^+||h5(i2-i2o)IU2(r)||li5i2o|U2(r). 

■"div ) -"div ) 



It remains to estimate ||*|| _ 1 . For <p e (F), let tpo denote the elementwise 

■^div 

average of (p. The Bramble-Hilbert Lemma yields 



||h H<p - ipo)\\L^{r) ^ Ivljji(p)' 
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which in conjunction with (28) implies 

(divr*T,iy?)i,r = y" divr^{(f - (po) dx ^ ||h5divr*||i2(r)|¥'|^i ^j,^ 

because of the norm equivalence for the discrete function Now, by definition 

(divr*,<^)ir i 
||divr*||^_i = sup — ^\\h-^Ih\\LHr), 



and 



Consequently 



(r) 



l*IL-i ^ ||h5Ho||i2(r) 



which together with (29) implies that 

Wh'-RoWmr) ^ ||u-C/|| 1 +||h5(il-i2o)||i2(r). 

■"div > 

Invoking the triangle inequality finally finishes the proof. □ 

It is important to realize the global nature of the above lower bound. This is due 
to the presence of integral operators Vk,Ak in the sesquilinear form a(-,-) which 
lead to a global norm for the error in (29) regardless of the support of 

In a very similar fashion, the following theorem can also be proven. 

_ 1 

Theorem 4.6 (Global lower bound for the curl residual). Let u G H^.^ (F) be the 

exact solution of (16) and U G i?To he its approximation defined by (17). Then, 
there exists a constant C3 > 0, only depending on shape regularity of T, such that 
the following bound holds 

C3||h5r|U2(r) < -1 ^ +||h^(r-ro)|U2(r). 

5. Conclusions 

In this paper we develop the first a posteriori error estimates for the electric 
field integral equation on polyhedra. We choose, for simplicity, to derive residual 
based error estimates but believe that our theory extends to other non-residual es- 
timators. We also choose to develop the theory for polyhedra, the most interesting 
and useful case in practice, but we expect the results to extend to smooth surfaces. 
For scattering problems on polyhedra, the solution u of the integral equation, or 
surface current, is not smooth whereas the regularity of the right-hand side / is 
dictated by the surface F because the incident wave is always smooth. This justi- 
fies our additional regularity assumption (21) which, coupled with the properties 
gradp(V/jdivrC/) G L^(F), curlr(Afet/) G i^^(F), allows us to evaluate the residuals 
R,r of (27) in i^(F) and thus avoid dealing with fractional Sobolev norms. We 
derive computable global upper and lower a posteriori bounds for the estimator (up 
to oscillation terms). In contrast to PDE, the estimator is global and due to the 
presence of the potentials V^, in the definition of the sesquilinear form. How- 
ever, the residuals R, r being evaluated in L'^{T) can be split elementwise and used 
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to drive an adaptive boundary element method (ABEM). The actual implementa- 
tion of ABEM for EFIE is rather delicate and is not part of the current discussion, 
which focusses on the derivation and properties of the estimators. 
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